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Abstract. Dyson- Schwinger equations furnish a Poincare covariant framework within 
which to study hadrons. A particular feature is the existence of a nonperturbative, sym- 
metry preserving truncation that enables the proof of exact results. Key to the DSE's 
efficacious application is their expression of the materially important momentum- 
dependent dressing of parton propagators at infrared length-scales, which is responsible 
for the magnitude of constituent-quark masses and the length-scale characterising con- 
finement in bound states. A unified quantitative description of light- and heavy-quark 
systems is achieved by capitalising on these features. 

1 Introduction 

This contribution provides an overview of one particular means by which a quan- 
titative and intuitive understanding of strong interaction phenomena can be 
attained. The broad framework is that of continuum strong QCD, by which 
I mean the continuum nonperturbative methods and models that can address 
these phenomena, especially those where a direct connection with QCD can be 
established, in one true limit or another. Naturally, everyone has a favourite tool 
and, in this connection, the Dyson-Schwinger equations (DSEs) are mine 1 . The 
framework is appropriate here because the last decade has seen a renaissance in 
its phenomenological application, with studies of phenomena as apparently un- 
connected as low-energy 7T7t scattering, B — > D* decays and the equation of state 
for a quark gluon plasma |2l,'il 1| . Indeed, the DSEs promise a single structure 
applicable to the gamut of strong interaction observables. 

Dyson-Schwinger equations provide a nonperturbative means of analysing a 
quantum field theory. Derived from a theory's Euclidean space generating func- 
tional, they are an enumerable infinity of coupled integral equations whose solu- 
tions are the n-point Schwinger functions (Euclidean Green functions), which are 
the same matrix elements estimated in numerical simulations of lattice-QCD. In 
theories with elementary fermions, the simplest of the DSEs is the gap equation, 
which is basic to studying dynamical symmetry breaking in systems as disparate 
as ferromagnets, superconductors and QCD. The gap equation is a good exam- 
ple because it is familiar and has all the properties that characterise each DSE. 
Its solution is a 2-point function (the fermion propagator) but its kernel involves 
higher rt-point functions; e.g., in a gauge theory, the kernel is constructed from 
the gauge-boson 2-point function and fcrmion-gauge-boson vertex, a 3-point 
function. In addition, while a weak-coupling expansion yields all the diagrams 
of perturbation theory, a self-consistent solution of the gap equation exhibits 
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nonperturbative effects unobtainable at any finite order in perturbation theory; 
e.g, dynamical chiral symmetry breaking (DCSB). 

The coupling between equations; namely, the fact that the equation for a 
given m-point function always involves at least one n > m-point function, neces- 
sitates a truncation of the tower of DSEs in order to define a tractable problem. 
It is unsurprising that the best known truncation scheme is just the weak cou- 
pling expansion which reproduces every diagram in perturbation theory. This 
scheme is systematic and valuable in the analysis of large momentum transfer 
processes because QCD is asymptotically free. However, it precludes the study of 
nonperturbative effects, and hence something else is needed for the investigation 
of strongly interacting systems and bound state phenomena. 

In spite of the need for a truncation, gap equations have long been used 
effectively in obtaining nonperturbative information about many-body systems 
as, e.g., in the Nambu-Gorkov formalism for superconductivity. The positive 
outcomes have been achieved through the simple expedient of employing the 
most rudimentary truncation, e.g., Hartree or Hartree-Fock, and comparing the 
results with observations. Of course, agreement under these circumstances is not 
an unambiguous indication that the contributions omitted are small nor that 
the model expressed in the truncation is sound. However, it does justify further 
study, and an accumulation of good results is grounds for a concerted attempt to 
substantiate a reinterpretation of the truncation as the first term in a systematic 
and reliable approximation. 

The modern application of DSEs, notably, comparisons with and predic- 
tions of experimental data, can properly be said to rest on model assumptions. 
However, those assumptions can be tested within the framework and also via 
comparison with lattice-QCD simulations, and the predictions are excellent. Fur- 
thermore, progress in understanding the intimate connection between symme- 
tries and truncation schemes has enabled exact results to be proved. Herein I 
will briefly explain recent phenomenological applications and the foundation of 
their success, and focus especially on the links the approach provides between 
light- and heavy-quark phenomena. It will become apparent that the momentum- 
dependent dressing of the propagators of QCD's elementary excitations is a 
fundamental and observable feature of strong QCD. 

The article is organised as follows: Sec.|21[p.E] - a review of DSE quiddities, 
especially in connection with the development of a nonperturbative, systematic 
and symmetry preserving truncation scheme, and the model-independent results 
whose proof its existence enables; Sec. [p. El ~ an illustration of the effica- 
cious application of DSE methods to light-meson systems and the connections 
that may be made with the results of lattice-QCD simulations; Sec. 0] [p. [33] - 
the natural extension of these methods to heavy-quark systems, with an expla- 
nation of the origin and derivation of heavy-quark symmetry limits and their 
confrontation with the real-world of finite quark masses; and Sec. [p. EE] - an 
epilogue. 
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2 Dyson-Schwinger Equations 
2.1 Gap Equation 

The simplest DSE is the gap equation, which describes how the propagation of 
a fermion is modified by its interactions with the medium being traversed. In 
QCD that equation assumes the form: 1 

S-^P) = Z 2 ((,A)i 1 - P + Z 4 (C, A) m(C) + S'ip, A), (1) 
wherein the dressed-quark self-energy is 

S'{p, A) = Z X (C, A)j\ 2 D^{p - q) S(q) r*(q,p). (2) 

Equations Q , @ constitute the renormalised DSE for the dressed-quark propa- 
gator. In Eq. (0), D llv {k) is the renormalised dressed-gluon propagator, r®(q;p) 
is the renormalised dressed-quark-gluon vertex and := j A d i q/(2Tr)' 1 repre- 
sents a translationally-invariant regularisation of the integral, with A the regu- 
larisation mass-scale. 2 In addition, Zi((,A), Z2(C,A) and Z4((,A) are, respec- 
tively, Lagrangian renormalisation constants for the quark-gluon vertex, quark 
wave function and quark mass-term, which depend on the renormalisation point, 
and the regularisation mass-scale, as does the gauge-independent mass renor- 
malisation constant, 

Z m (C 2 , A 2 ) = Z A (C 2 ,A 2 ) Z^(C 2 ,A 2 ) , (3) 
whereby the renormalised running-mass is related to the bare mass: 

m(0 = Z- n 1 (( 2 ,A 2 )m hm (A)- (4) 

When £ is very large the running-mass can be evaluated in perturbation theory, 
which gives 

jn 

m(CH (lnT7^)-' 7m = 12/(33 ~ 2Nf) - ^ 

Here Nf is the number of current-quark flavours that contribute actively to the 
running coupling, and Aqcb and m are renormalisation group invariants. 

The solution of Eq. is the dressed-quark propagator and takes the form 

S-\p) = i~i - P A(p 2 , C 2 ) + B(p 2 , C 2 ) ee 1 [ ll -p + M(p 2 )] . (6) 



1 A Euclidean metric is employed throughout, wherewith the scalar product of two 
four vectors is a ■ b — J^- = i a i^i an< A 1 employ Hermitian Dirac-7 matrices that obey 
{7n,7f} = 2(5 M!/ and tr 757^7,, 7^7,7 = — 4e M „ pCT , £1234 = 1. 

2 Only with a translationally invariant regularisation scheme can Ward-Takahashi 
identities be preserved, something that is crucial to ensuring vector and axial- vector 
current conservation. The final stage of any calculation is to take the limit A — > 00. 
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It is obtained by solving the gap equation subject to the renormalisation condi- 
tion that at some large spacelike £ 2 

S-'ip^^^H-p + miQ. (7) 

The observations made in in the Introduction are now manifest. The gap 
equation is a nonlinear integral equation for S (p) and can therefore yield much- 
needed nonperturbative information. However, the kernel involves the two-point 
function D^ v (k) and the three-point function r®(q;p). The equation is conse- 
quently coupled to the DSEs these functions satisfy and hence a manageable 
problem is obtained only once a truncation scheme is specified. 



2.2 Nonperturbative Truncation 

To understand why Eq. Q is called a gap equation, consider the chiral limit, 
which is readily defined |S] because QCD exhibits asymptotic freedom and im- 
plemented in the gap equation by employing 6 

Z 2 (( 2 ,A 2 )m hnl (A) =0, yl»C. (8) 

It is noteworthy that for finite ( and A — > oo, the left hand side (l.h.s.) of Eq. |(HJ 
is identically zero, by definition, because the mass term in QCD's Lagrangian 
density is renormalisation-point-independent. The condition specified in Eq. ©, 
on the other hand, effects the result that at the (perturbative) renormalisation 
point there is no mass-scale associated with explicit chiral symmetry breaking, 
which is the essence of the chiral limit. An equivalent statement is that one 
obtains the chiral limit when the renormalisation-point-invariant current-quark 
mass vanishes; namely, rh = in Eq. JSJ). In this case the theory is chirally 
symmetric, and a perturbative evaluation of the dressed-quark propagator from 
Eq. (P) gives 




viz., the perturbative mass function is identically zero in the chiral limit. It 
follows that there is no gap between the top level in the quark's filled negative- 
energy Dirac sea and the lowest positive energy level. 

However, suppose one has at hand a truncation scheme other than pertur- 
bation theory and that subject to this scheme Eq. JQ) possessed a chiral limit 
solution B°(p 2 ) 0. Then interactions between the quark and the virtual quanta 
populating the ground state would have nonperturbatively generated a mass gap. 
The appearance of such a gap breaks the theory's chiral symmetry. This shows 
that the gap equation can be an important tool for studying DCSB, and it has 
long been used to explore this phenomenon in both QED and QCD £Q. 

The gap equation's kernel is formed from a product of the dressed-gluon 
propagator and dressed-quark-gluon vertex but in proposing and developing a 
truncation scheme it is insufficient to focus only on this kernel |b!7| . The gap 
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equation can only be a useful tool for studying DCSB if the truncation itself 
does not destroy chiral symmetry. 

Chiral symmetry is expressed via the axial-vector Ward-Takahashi identity: 

P ti r^(k;P) = S- 1 {k + )i l5 +i l5 S- 1 (k_), k ± =k±P/2, (10) 

wherein r^^{k\P) is the dressed axial- vector vertex. This three-point function 
satisfies an inhomogeneous Bethe-Salpeter equation (BSE): 

[r 5 p(k; P)} tu = Z 2 [ 757/1 L + [ A [S(q+)r 5 ,(q; P)S(q-)]„K£(q, k; P) , (11) 

Jq 

in which K(q, k] P) is the fully-amputated quark-antiquark scattering kernel, and 
the colour-, Dirac- and flavour-matrix structure of the elements in the equation 
is denoted by the indices r, s, t, u. The Ward-Takahashi identity, Eq. (|1C)|) . entails 
that an intimate relation exists between the kernel in the gap equation and that 
in the BSE. (This is another example of the coupling between DSEs.) There- 
fore an understanding of chiral symmetry and its dynamical breaking can only 
be obtained with a nonperturbative truncation scheme that preserves this rela- 
tion, and hence guarantees Eq. (|10|) without a fine-tuning of model-dependent 
parameters. 



Rainbow-ladder truncation. At least one such scheme exists Its leading- 
order term is the so-called renormalisation-group-improved rainbow- ladder trun- 
cation, whose analogue in the many body problem is an Hartree-Fock truncation 
of the one-body (Dyson) equation combined with a consistent ladder-truncation 
of the related two-body (Bethe-Salpeter) equation. To understand the origin of 
this leading-order term, observe that the dressed-ladder truncation of the quark- 
antiquark scattering kernel is expressed in Eq. I|ll|l via 

[L(q,k;P)]t?' [r 5fi (q;P)] u , t , := [S(q+)r 5fi (q; P)S(q„)] sr K%(q, k; P) 
= -g 2 (C 2 )D pa (k-q) 

x [r^k + ,q + )S(q + )] tt , [S(q.)r^,k.)] u , u [r 5fi (q; P)] t , u , (12) 

wherein I have only made explicit the renormalisation point dependence of the 
coupling. One can exploit multiplicative renormalisability and asymptotic free- 
dom to demonstrate that on the kinematic domain for which Q 2 := (k — q) 2 ~ 
k 2 ~ q 2 is large and spacelike 



[L(q,k;P)]%' =-4™(Q 2 )D 1 ™(Q) 



y7pS frc °(^ 



S ficc (q-) y7. 



(13) 

where a(Q 2 ) is the strong running coupling and, e.g., S ° is the free quark prop- 
agator. It follows that on this domain the r.h.s. of Eq. describes the lead- 
ing contribution to the complete quark-antiquark scattering kernel, K[£(q, k; P), 
with all other contributions suppressed by at least one additional power of l/Q 2 . 



6 Craig D. Roberts 



The renormalisation-group-improved ladder-truncation supposes that 



K™(q, k- P) = -47ra(Q 2 ) D^(Q) 



"A a 






Y 7p_ 


ts 





(14) 



is also a good approximation on the infrared domain and is thus an assumption 
about the long-range (Q 2 < 1 GeV 2 ) behaviour of the interaction. Combining 
Eq. Q14JI with the requirement that Eq. Ijl()|l be automatically satisfied leads to 
the renormalisation-group-improved rainbow-truncation of the gap equation: 



>' V> = Z-i(h + / 4tt a(Q 2 ) D^(p - q)^ S(q) y 7 * . (15) 



This rainbow-ladder truncation provides the foundation for an explanation of a 
wide range of hadronic phenomena 0] . 



2.3 Systematic Procedure 

The truncation scheme of Ref. [H] is a dressed-loop expansion of the dressed- 
quark-gluon vertices that appear in the half-amputated dressed-quark-antiquark 
scattering matrix: S 2 K, a renormalisation-group invariant [§]. All n-point func- 
tions involved thereafter in connecting two particular quark-gluon vertices are 
fully dressed. The effect of this truncation in the gap equation, Eq. is re- 
alised through the following representation of the dressed-quark-gluon vertex, 

ir° = ±\T lt = iT lt : 

Z 1 r ll {k,p) = 7m + ^- J g 2 D p<J { P - £) lp S(£ + k-p)^S(£) 7<T 
N f A 

+ Ti 9 D «'*( i ) D T'T(e+k-p)jr'S( P -£) 7l ,,r% fl (£,-k,k-p) + [...]. 

(16) 

Here r 3g is the dressed-three-gluon vertex and it is readily apparent that the 
lowest order contribution to each term written explicitly is 0(<? 2 ). The ellipsis 
represents terms whose leading contribution is 0(g 4 ); viz., the crossed-box and 
two-rung dressed-gluon ladder diagrams, and also terms of higher leading-order. 

This expansion of S 2 K, with its implications for other n-point functions, 
yields an ordered truncation of the DSEs that guarantees, term-by-term, the 
preservation of vector and axial- vector Ward-Takahashi identities, a feature that 
has been exploited [5I1UI11) to establish exact results in QCD. It is readily seen 
that inserting Eq. (|16f> into Eq. £TJ provides the rule by which the rainbow-ladder 
truncation can be systematically improved. 

Planar vertex. The effect of the complete vertex in Eq. Ijl6|) on the solutions of 
the gap equation is unknown. However, insights have been drawn from a study 
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ri(k,p) = 



Fig. 1. Integral equation for a planar dressed-quark-gluon vertex obtained by neglect- 
ing contributions associated with explicit gluon self-interactions. Solid circles indicate 
fully dressed propagators. The vertices are not dressed. (Adapted from Ref. 0.) 

9 of a more modest problem obtained by retaining only the sum of dressed- 
gluon ladders; i.e., the vertex depicted in Fig.Q] The elucidation is particularly 
transparent when one employs |12| 

2V(fc) := g 2 D^(k) = (s^ - (2vr) 4 Q 2 <5 4 (fc) (17) 

for the dressed-gluon line, which defines an ultraviolet finite model so that the 
regularisation mass-scale can be removed to infinity and the renormalisation 
constants set equal to one. 3 This model has many positive features in com- 
mon with the class of renormalisation-group-improved rainbow-ladder models 
and its particular momentum-dependence works to advantage in reducing in- 
tegral equations to algebraic equations with similar qualitative features. There 
is naturally a drawback: the simple momentum dependence also leads to some 
model-dependent artefacts, but they are easily identified and hence not cause 
for concern. 

The general form of the dressed-quark gluon vertex involves twelve distinct 
scalar form factors but using Eq. (jT7|) only I),(p) := r^{p,p) contributes to 
the gap equation, which considerably simplifies the analysis. The summation 
depicted in Fig. ^ is expressed via 

r» (?) = i» + li P s( P ) r M ( P ) s( P ) 7p , (is) 

which supports a solution 

r^(p) ="i(p 2 )7m+ ^ 2 (p 2 )l ■ PPf,- otaiP^iPn- (19) 
One can re-express this vertex as 

oo oo 

r M (p) = J2 r » = E K(P 2 ) ^ + 4b 2 ) 7 • PP, - 4(p 2 ) iP/J > (20) 

Z=0 2=0 

3 The constant Q sets the model's mass-scale and using Q = 1 simply means that all 
mass-dimensioned quantities are measured in units of Q. 
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where the superscript enumerates the order of the iterate: r^ =0 is the bare vertex, 

al = 1 , a° 2 = = a% ; (21) 

-T^ =1 is the result of inserting this into the r.h.s. of Eq H18|l to obtain the one- rung 
dressed-gluon correction; i~^ =2 is the result of inserting i~^ =1 , and is therefore 
the two-rung dressed-gluon correction; etc. A key observation is that each 
iterate is related to its precursor via a simple recursion relation and, substituting 
Eq. (J20J, that recursion yields (s = p 2 ) 

(a\ +1 (s)\ 

a> +1 (s):= ai +1 (s)\=0(s;A,B) a i (s), (22) 

WW 
1 1 f- A 

0(s; A, B) = - 2A 2 sA 2 - B 2 2AB \ ( ->:! ! 

\AAB AsAB 2(B 2 -sA 

A = sA 2 (s) + B 2 (s). It follows that 

«=ff:^a = T 4^« (24) 




and hence, using Eq. 

AA 

ai 



1 + AA 



-8A 2 1 + 2A 

2 1 + 2{B 2 -sA 2 )~8A 2 1 + 4, A ' 1 ' 

-8AB 

" 3 ~ l + 2(£ 2 -sv4 2 )-8Z\ 2 ' 

The recursion relation thus leads to a closed form for the gluon-ladder-dressed 
quark-gluon vertex in Fig.^ viz., Eqs. I|19[) . I|25(l . Its momentum-dependence is 
determined by that of the dressed-quark propagator, which is obtained by solving 
the gap equation, itself constructed with this vertex. Using Eq. Ijl7|l . that gap 
equation is 

S- 1 (p)=i>y-p + m + 'y M S(p)r ll (p) (26) 
whereupon the substitution of Eq. (|19|l gives 

M») = 1 + sA2 \ B2 [A {2ax - sa 2 ) -Ba 3 ], (27) 

B ( s ) = m + a 2 \ n 2 i B ( 4 «i + s «2) - sA a 3 ] . (28) 
sA z + B z 

Equations l|27|l . I|28l) . completed using Eqs. (|25l) . form a closed algebraic system. 
It can easily be solved numerically, and that yields simultaneously the complete 
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Fig. 2. A(a), M(s) obtained from Eqs. ^J, g3, with m = 0.023 (solid line). 
All dimensioned quantities are expressed in units of Q in Eq. $17\ . For comparison, 
the results obtained with the zeroth-order vertex (dotted line) and the one-loop vertex 
(dashed line) are also plotted. (Adapted from Ref. 

gluon-ladder-dressed vertex and the propagator for a quark fully dressed via 
gluons coupling through this nonperturbative vertex. Furthermore, it is apparent 
that in the chiral limit, m — 0, a realisation of chiral symmetry in the Wigner- 
Weyl mode, which is expressed via the B = solution of the gap equation, is 
always admissible. This is the solution anticipated in Eq. JJjJ. 

The chiral limit gap equation also admits a Nambu-Goldstone mode solution 
whose p 2 ~ properties are unambiguously related to those of the m ^ so- 
lution, a feature also evident in QCD A complete solution of Eq. I)26|l is 
available numerically, and results for the dressed-quark propagator are depicted 
in Fig. [21 It is readily seen that the complete resummation of dressed-gluon lad- 
ders gives a dressed-quark propagator that is little different from that obtained 
with the one-loop-corrected vertex; and there is no material difference from the 
result obtained using the zeroth-order vertex. Similar observations apply to the 
vertex itself. The scale of these modest effects can be quantified by a comparison 
between the values of M(s = 0) = B(0)/A(0) calculated using vertices dressed 
at different orders: 





N = 


N = 1 


N = 2 


N = oo 


M(0) 


1 


1.105 


1.115 


1.117 



(29) 



The rainbow truncation of the gap equation is accurate to within 12% and adding 
just one gluon ladder gives 1% accuracy. It is important to couple this with an 
understanding of how the vertex resummation affects the Bethe-Salpeter kernel. 

Vertex-consistent Bethe-Salpeter kernel. The renormalised homogeneous 
BSE for the quark-antiquark channel denoted by M can be expressed 

[r M (k- P)] tu = I [ XM {q; P)] sr [K(k, q; P)]£ , (30) 

Jq 
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Fig. 3. Homogeneous BSE, Eq. l-'ilH . Filled circles: dressed propagators or vertices; 
K is the dressed-quark-antiquark scattering kernel. A systematic truncation of S 2 K is 
the key to preserving Ward-Takahashi identities {SH5] . (Adapted from Ref. 



where: Pm (fc; P) is the meson's Bethe-Salpeter amplitude, k is the relative mo- 
mentum of the quark-antiquark pair, P is their total momentum; and 

X M (k; P) = S(k+) r M (k; P) . (31) 

Equation (|3Ufl. depicted in Fig. [31 describes the residue at a pole in the solution 
of an inhomogeneous BSE; e.g., the lowest mass pole solution of Eq. I|ll|) is 
identified with the pion. 4 

I noted on p. that the automatic preservation of Ward-Takahashi identities 
in those channels related to strong interaction observables requires a conspiracy 
between the dressed-quark-gluon vertex and the Bethe-Salpeter kernel |8I14| . A 
systematic procedure for building that kernel follows [5] from the observation 
[TU that the gap equation can be expressed via 

™ = o, 

where r[S] is a Cornwall- Jackiw-Tomboulis- like effective action |TS]. The Bethe- 
Salpeter kernel is then obtained via an additional functional derivative: 

KZ = - 5 -^. (33) 

With the vertex depicted in Fig. ^ the n-th order contribution to the kernel 
is obtained from the n-loop contribution to the self energy: 

£ n (p) = ~ [ A V^(p-q)l a ltl S(q)l a r?(q,p). (34) 
Jq 

Since P M (p, q) itself depends on S then Eq. (fT?T?f> yields the Bethe-Salpeter kernel 
as a sum of two terms and hence Eq. (f3T)|) assumes the form 



r M (k;P)= [ V„ u (k-q)l a j„ 

J a 



XM (q; P) l a r v (q-,k-) + S(q+) A a M „(q, k; P) , 

(35) 

4 The canonical normalisation of a Bethe-Salpeter amplitude is fixed by requiring that 
the bound state contribute with unit residue to the fully-amputated quark-antiquark 
scattering amplitude: M = K + K(SS)K + [. . .]. See, e.g., Ref. [TS]. 
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Fig. 4. Upper panel: BSE, Eq. (I35L which is valid whenever 7^ can be obtained via a 
recursion relation. Lower panel: Recursion relation for A 1 ^. (Adapted from Ref. [5].) 

where I have used the mnemonic 



Equation l(35jl is depicted in the upper panel of Fig. The first term is 
instantly available once one has an explicit form for J 1 ™ and the second term, 
identified by the shaded box in Fig.0] can be obtained [5] via the inhomogeneous 
recursion relation depicted in the figure's lower panel. Combining these figures, it 
is apparent that to form the Bethe-Salpeter kernel the free gluon line is attached 
to the upper dressed-quark line. Consequently, the first term on the r.h.s. of 
the lower panel in Fig. ^invariably generates crossed gluon lines; viz., nonplanar 
contributions to the kernel. The character of the vertex-consistent Bethe-Salpeter 
kernel is now clear: it consists of countably many contributions, a subclass of 
which are crossed-ladder diagrams and hence nonplanar. Only the rainbow gap 
equation, obtained with i = in Eq. ((2()|) . yields a planar vertex-consistent 
Bethe-Salpeter kernel, namely the ladder kernel of Eq. 114(1 . In this case alone 
is the number of diagrams in the dressed-vertex and kernel identical. Otherwise 
there are always more terms in the kernel. 

Solutions for the n- and p-mesons. I have recapitulated on a general proce- 
dure that provides the vertex-consistent channel-projected Bethe-Salpeter ker- 
nel once r™ and the propagator functions; A, B, are known. That kernel must 
be constructed independently for each channel because, e.g., A a Mv depends on 
Xm(<?; P)- As with the study of the vertex, an elucidation of the resulting BSEs' 
features is simplified by using the model of Eq. (fT7|l . for then the Bethe-Salpeter 
kernels are finite matrices [cf. (1 — O) -1 in Eq. 1|24|) ] and the homogeneous BSEs 
are merely linear, coupled algebraic equations. 

Reference |3j describes in detail the solution of the coupled gap and Bethe- 
Salpeter equations for the tt- and p-mesons. Herein I focus on the results, which 
are summarised in Tabled It is evident that, irrespective of the order of the 
truncation; viz., the number of dressed gluon rungs in the quark-gluon vertex, the 



oo 




(36) 
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Table 1. Calculated n and p meson masses, in GeV, quoted with Q = 0.48 GeV, in 
which case m — 0.023 5 = HMeV. n is the number of dressed-gluon rungs retained 
in the planar vertex, see Fig. Q and hence the order of the vertex-consistent Bethe- 
Salpeter kernel: the rapid convergence of the kernel is apparent from the tabulated 
results. (Adapted from Ref. 0.) 
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pion is massless in the chiral limit. This is in spite of the fact that it is composed 
of heavy dressed-quarks, as is clear in the calculated scale of the dynamically 
generated dressed-quark mass function: see Fig. |3 M(0) w Q « 0.5 GeV. These 
observations emphasise that the masslessness of the 7r is a model-independent 
consequence of consistency between the Bethe-Salpeter kernel and the kernel 
in the gap equation. Furthermore, the bulk of the p-n mass splitting is present 
for m — and with the simplest (n = 0; i.e., rainbow-ladder) kernel, which 
demonstrates that this mass difference is driven by the DCSB mechanism. It 
is not the result of a carefully contrived chromo-hyperfine interaction. Finally, 
the quantitative effect of improving on the rainbow-ladder truncation; namely, 
including more dressed-gluon rungs in the gap equation's kernel and consistently 
improving the kernel in the Bethe-Salpeter equation, is a 10% correction to the 
vector meson mass. Simply including the first correction (viz., retaining the first 
two diagrams in Fig. ^) yields a vector meson mass that differs from the fully 
resummed result by rs 1%. The rainbow-ladder truncation is clearly accurate in 
these channels. 

Comments. While I have described results obtained with a rudimentary in- 
teraction model in order to make the construction transparent, the procedure 
is completely general. However, the algebraic simplicity of the analysis is nat- 
urally peculiar to the model. With a more realistic interaction, the gap and 
vertex equations yield a system of twelve coupled integral equations. The Bethe- 
Salpeter kernel for any given channel then follows as the solution of a determined 
integral equation. 

The material reviewed covers those points in the construction of Refs. |8I9| 
that bear upon the fidelity of the rainbow-ladder truncation in pairing the gap 
equation and Bethe-Salpeter equations for the vector and flavour non-singlet 
pseudoscalar mesons. The error is small. In modelling it is therefore justified to 
fit one's parameters to physical observables at this level in these channels and 
then make predictions for other phenomena involving vector and pseudoscalar 
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bound states in the expectation they will be reliable. That approach has been 
successful, as illustrated in Ref. 0|. 

Lastly, the placement of the rainbow-ladder truncation as the first term in a 
procedure that can methodically be improved explains why this truncation has 
been successful, the boundaries of its success, why it has failed outside these 
boundaries, and why sorting out the failures won't undermine the successes. 

2.4 Selected Model-Independent Results 

In the hadron spectrum the pion is identified as both a Goldstone mode, associ- 
ated with DCSB, and a bound state composed of constituent u- and ci-quarks, 
whose effective mass is ~ 300 — 500 MeV. Naturally, in quantum mechanics, one 
can fabricate a mass operator that yields a bound state whose mass is much 
less than the sum of the constituents' masses. However, that requires fine tun- 
ing and, without additional fine tuning, such models predict properties for spin- 
and/or isospin-flip relatives of the pion which conflict with experiment. A cor- 
rect resolution of this apparent dichotomy is one of the fundamental challenges 
to establishing QCD as the theory underlying strong interaction physics, and 
the DSEs provide an ideal framework within which to achieve that end, as I 
now explain following the proof of Ref. |S] . It cannot be emphasised too strongly 
that the legitimate understanding of pion observables; including its mass, decay 
constant and form factors, requires an approach to contain a well-defined and 
valid chiral limit. 



Proof of Goldstone's Theorem. Consider the BSE expressed for the isovector 
pseudoscalar channel: 



lxUq;P)}srK™(q,k;P), 



(37) 



with xi(l-P) = S(q + )ri(q;P)S(q^) obvious from Eq. (J3TJ) and j labelling 
isospin, of which the solution has the general form 



ri(k;P)=T> 75 



iE v (k;P)+i-PF v (k;P) 



7 -kk-PG„(k;P) + a flv k^P„ H^(k; P) 



(38) 



It is apparent that the dressed-quark propagator, the solution of Eq. (Q, is an 
important part of the BSE's kernel. 

Chiral symmetry and its dynamical breaking are expressed in the axial- vector 
Ward-Takahashi identity, Eq. (|10|l . which involves the axial- vector vertex: 



r£jk-,p) 



757m - 



+ / H^P^srKZiq^P), (39) 
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that has the general form 

%(k; P) = y 75 [^F R (k; P) + 7 • kk M G R {k; P) - k v H R (k; P)} 

+ n,(k;P) + 1 ^- I ^(k;P), (40) 

where F R , G R , H R and % are regular as P 2 -» -m 2 , P„f^(k;P) ~ 0(P 2 ) 
and (j>> (k; P) has the structure depicted in Eq. (|38|) . Equation (|40|l admits the 
possibility of at least one pole term in the vertex but does not require it. 

Substituting Eq. 140(1 into l|39|) and equating putative pole terms, it is clear 
that, if present, (j^{k;P) satisfies Eq. (|37fl . Since this is an eigenvalue problem 
that only admits a solution for P 2 — —m 2 , it follows that <j>>(k\P) 

is nonzero solely for P 2 = —m 2 and the pole mass is m 2 ^ = m 2 . Hence, if 
K supports such a bound state, the axial-vector vertex contains a pion-pole 
contribution. Its residue, r A , however, is not fixed by these arguments. Thus Eq. 
(|40() becomes 

r^(k- P) = y 75 [^F R (k; P) + 7 • kk M G R {k; P) - k v H R (k; P)} 

+n^p)+^^m-,p). (4i) 

Consider now the chiral limit axial- vector Ward-Takahashi identity Eq. (|10|) . 
If one assumes m 2 = in Eq. Q4ip. substitutes it into the l.h.s. of Eq. I|10|) along 
with Eq. © on the right, and equates terms of order (P„)° and P„, one obtains 
the chiral-limit relations [5] 

r A E 7r (k;0) = B(k 2 ), F R (k; 0) + 2 r A F„{k; 0) = A(k 2 ) , 
G R (k;0) + 2r A G 7r (k;0)=2A'(k 2 ), H R (k; 0) + 2r A H T (k; 0) = . 1 } 

I have already explained that B(k 2 ) = in the chiral limit [remember Eq. 
and that a B(k 2 ) ^ solution of Eq. (|TJ in the chiral limit signals DCSB. Indeed, 
in this case QB] 

M(p2) 2^ (~ ^ , (43 ) 

p 2 U\n p 2 /A 2 QCB 



where (qq}° is the renormalisation-point-independent vacuum quark condensate 
|17|. Furthermore, there is at least one nonperturbative DSE truncation scheme 
that preserves the axial-vector Ward-Takahashi identity order by order. Hence 
Eqs. (|42|l are exact quark-level Goldberger-Treiman relations, which state that 
when chiral symmetry is dynamically broken: 

(i) . the homogeneous isovector pseudoscalar BSE has a massless solution; 

(ii) . the Bethe-Salpeter amplitude for the massless bound state has a 

term proportional to 75 alone, with 0) completely determined 

by the scalar part of the quark self energy in addition to other 
pseudoscalar Dirac structures, F n , G^ and H^, that are nonzero; 

(iii) . and the axial- vector vertex is dominated by the pion pole for P 2 ~ 0. 
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The converse is also true. Hence DCSB is a sufficient and necessary condition 
for the appearance of a massless pseudoscalar bound state (of what can be very- 
massive constituents) that dominates the axial- vector vertex for P 2 as 0. 



Mass Formula. When chiral symmetry is explicitly broken the axial-vector 
Ward-Takahashi identity becomes: 

JVt(*i P) = 5- 1 (fc + ) l75 y + ?7 5y 5-!(fc_) - 2* m(C) P|(fc; P), (44) 
where the pseudoscalar vertex is obtained from 

ri(k-,p) 

As argued in connection with Eq. (|39|) , the solution of Eq. I|45|l has the form 

iri (k; P) = y 75 [iE?{k- P)+ 1 -PF?+ 1 -kk-PG p R {k-P) 

+ k^P v ffg (fc; P)] + ? ' P 2 Pi (fc; P) , (46) 

"+" TO ir 

where P^, P^, G R and Pj^ are regular as P 2 — > — m 2 ; i.e., the isovector pseu- 
doscalar vertex also receives a contribution from the pion pole. In this case 
equating pole terms in the Ward-Takahashi identity Eq. H44J) . entails [H] 

r A ml = 2m(0MC). ( 47 ) 

This, too, is an exact relation in QCD. Now it is important to determine the 
residues ta and rp. 

Study of the renormalised axial- vector vacuum polarisation shows jS] : 

ta 6V P, = U S" P, = Z 2 tr / ir i 7 57 AI 5(g + )Pi( (Z ; P)S(<?_) , (48) 

Jq 

where the trace is over colour, Dirac and flavour indices; i.e., the residue of the 
pion pole in the axial- vector vertex is the pion decay constant. The factor of 
Z2 on the r.h.s. in Eq. I|48|l is crucial: it ensures the result is gauge invariant, 
and cutoff and renormalisation-point independent. Equation (|48J) is the exact 
expression in quantum field theory for the pseudovector projection of the pion's 
wave function on the origin in configuration space. 

A close inspection of Eq. (|45|l . following its re-expression in terms of the 
renormalised, fully-amputated quark-antiquark scattering amplitude: M — K + 
K(SS)K + . . ., yields 

ir P 5« = Z A tr j A lr l l5 S(q+)rUq; P)S(q-) , (49) 

Jq 



75" 



+ 



K£(q,k',P). 



(45) 
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wherein the dependence of Z<± on the gauge parameter, the regularisation mass- 
scale and the renormalisation point is exactly that required to ensure: 1) rp is 
finite in the limit A — > oo; 2) rp is gauge-parameter independent; and 3) the 
renormalisation point dependence of rp is just such as to guarantee the r.h.s. 
of Eq. H47|) is renormalisation point independent. Equation l)49|) expresses the 
pseudoscalar projection of the pion's wave function on the origin in configuration 
space. 

Focus for a moment on the chiral limit behaviour of Eq. I|49(l whereat, using 
Eqs. pSjl. JUJl, one finds readily 

- (qq)° C = Ur° P (0 = Z 4 (C, A) N c tr D / S fn=0 (q) . (50) 

Jq 

Equation (|50|l is unique as the expression for the chiral limit vacuum quark 
condensate. 5 It is ^-dependent but independent of the gauge parameter and the 
regularisation mass-scale, and Eq. 1)50(1 thus proves that the chiral-limit residue 
of the pion pole in the pseudoscalar vertex is (— (qq)®) / ■ Now Eqs. 1)47(1 . 150JI 
yield 

(f^ml =-2m(C)y» + 0(m 2 ), (51) 

where f® is the chiral limit value from Eq. I(48|l . Hence what is commonly known 
as the Gell-Mann-Oakes-Renner relation is a corollary of Eq. I|47|l . 

One can now understand the results in Table ^ a massless bound state of 
massive constituents is a necessary consequence of DCSB and will emerge in any 
few-body approach to QCD that employs a systematic truncation scheme which 
preserves the Ward-Takahashi identities. 

Upon review it will be apparent that Eqs. (|47Jl - (|49|l are valid for any 
values of the current-quark masses, and the generalisation to Nf quark flavours 
is 

flrnl = -{qq)^M%, (52) 

= m^ qi + toJ> 2 is the sum of the current-quark masses of the meson's con- 
stituents; 

f A 

f H P^ = Z 2 tr \{T H ) T lbl ^S{q + )r H {q-P)S{q„), (53) 
Jq 

with S — diag(S' u , Sd, S s , . . .), T H a flavour matrix specifying the meson's quark 
content, e.g., T n = ^(X 1 + i\ 2 ), {A 4 } are iV/-flavour generalisations of the 
Gell-Mann matrices; and 

f A 

(qq)f = -f H ri = if H Z A tx ±(T H ) T l5 S(q + ) r H (q; P) 5(g_) . (54) 

Jq 

NB. Equation (|50|l means that in the chiral limit (qq)^ — > (qq)® and hence (qq)^ 
has been called an in-hadron condensate. 



5 The trace of the massive dressed-quark propagator is not renormalisable and hence 
there is no unique definition of a massive-quark condensate |17|. 
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The formulae reviewed in this Section also yield model-independent corollar- 
ies for systems involving heavy-quarks, as I relate in Sec. 



3 Basis for a Description of Mesons 

The renormalisation-group-improved rainbow-ladder truncation has long been 
employed to study light mesons, and in Sees. 12.21 12. 31 it was shown to be a quan- 
titatively reliable tool for vector and flavour nonsinglet pseudoscalar mesons. In 
connection with Eqs. I|14l) . (|15ll . I argued that the truncation preserves the ultra- 
violet behaviour of the quark-antiquark scattering kernel in QCD but requires an 
assumption about that kernel in the infrared; viz., on the domain Q 2 < 1 GeV 2 , 
which corresponds to length-scales > 0.2 fm. The calculation of this behaviour 
is a primary challenge in contemporary hadron physics and there is progress 
|7I18I19I20I21I22I215| . However, at present the most efficacious approach is to 
model the kernel in the infrared, which enables quantitative comparisons with 
experiments that can be used to inform theoretical analyses. The most exten- 
sively applied model is specified by [21] 



In 



2 

QCD 



F(Q 2 ) , (55) 



in Eqs. (O, G3. Here, T{Q 2 ) = [1 - exp(-Q 2 /[4m 2 ])]/Q 2 , m t = 0.5 GeV; 
t = e 2 - 1; 7m = 12/25; and [25] A QC o = a!£L = 0.234 GeV. This simple form 
expresses the interaction strength as a sum of two terms: the second ensures that 
perturbative behaviour is preserved at short-range; and the first makes provision 
for the possibility of enhancement at long-range. The true parameters in Eq. 155JI 
are D and u>, which together determine the integrated infrared strength of the 
rainbow-ladder kernel; i.e., the so-called interaction tension, a A However, 
I emphasise that they are not independent: in fitting to a selection of observ- 
ables, a change in one is compensated by altering the other; e.g., on the domain 
ijj G [0.3, 0.5] GeV, the fitted observables are approximately constant along the 
trajectory [7] 

ujD = (0.72 GeV) 3 . (56) 

Hence Eq. (|55J) is a one-parameter model. This correlation: a reduction in D 
compensating an increase in cj, ensures a fixed value of the interaction tension. 



3.1 Rainbow Gap Equation 

Equations l|15fl and (|55|l provide a model for QCD's gap equation and in hadron 
physics applications one is naturally interested in the nonperturbative DCSB 
solution. A familiar property of gap equations is that they only support such a 
solution if the interaction tension exceeds some critical value. In the present case 
that value is erf ~ 2.5 GeV/fm JSj- This amount of infrared strength is sufficient 
to generate a nonzero vacuum quark condensate but only just. An acceptable 
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Table 2. Comparison of experimental values with results for ir and K observables 
calculated using the rainbow-ladder interaction specified by Eq. I|55|l . quoted in MeV. 
The model's sole parameter and the current-quark masses were varied to obtain these 
results, and the best fit parameter values are given in Eqs. 1571 1 . 158H . Predictions 
for analogous vector meson observables are also tabulated. NB. A charged particle 
normalisation is used for in Eq. 1611 . which differs from that in Eq. 1531 by a 
multiplicative factor of y/2. (Adapted from Ref. |24|.1 
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description of hadrons requires a A ~ 25 GeV/fm |Hj and that is obtained with 

El 

D = (0.96 GeV) 2 . (57) 

This value of the model's infrared mass-scale parameter and the two current- 
quark masses 

m u (lGeV) = 5.5 MeV, m s (l GeV) = 125 MeV , (58) 
defined using the one-loop expression 

to evolve m u (19GeV) = 3.7 MeV and TO s (19GeV) = 85 MeV, were obtained in 
Ref. '21] by requiring a least-squares fit to the tt- and if-meson observables listed 
in Table [21 The procedure was straightforward: the rainbow gap equation [Eqs. 
0, JT5J, {55)] was solved with a given parameter set and the output used to 
complete the kernels in the homogeneous ladder BSEs for the n- and if-mesons 
[Eqs. O, JT5J, ® with T j for the tt channel and r-? -> T K+ = i(A 4 +iA 5 ) 
for the if]. These BSEs were solved to obtain the tt- and if-meson masses, and 
the Bethe-Salpeter amplitudes. Combining this information delivers the leptonic 
decay constants via Eq. 1)53(1 . This was repeated as necessary to arrive at the 
results in Table The model gives a vacuum quark condensate 

-<OT>?GeV= (0-242 GeV) 3 , (60) 

calculated from Eq. (|50|l and evolved using the one-loop expression in Eq. (|59(l . 

With the model's single parameter fixed, and the dressed-quark propagator 
determined, it is straightforward to compose and solve the homogeneous BSE 
for vector mesons. This yields predictions, also listed in Table for the vector 
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— DSE soln, ( 55) and ( 57) 
DSE soln, ( 55) and ( 62) 
• lattice Z(p) 
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Fig. 5. Left panel - wave function renormalisation: solution of the gap equation using 
Eqs. 1551 . 157H (solid curve); solution using Eqs. 15511 . 1621 (dashed- curve); quenched 
lattice-QCD simulations obtained with m = 0.036/a ~ 60MeV (data). The DSE 
study used a renormalisation point £ = 19GeV and a current-quark mass 0.6 ml GcV 
[Eq. 1581 1. to enable a direct comparison with the lattice data. Right panel - mass 
function. (Adapted from Ref. |7j.) 



meson masses and electroweak decay constants ^T] 

flMl = \Z 2 te { A (T H ) T ^S{ q+ )r»{q;P)S(q-), (61) 

Jq 

where Mj£ is the meson's mass and P^r^(q;P) = for P 2 = —(My) 2 ; i.e., 
the Bethe-Salpeter amplitude is transverse, characterises decays such as p — > 
e+e~, r — * K*v T . 

Given the discussion in Sec. [3 the phcnomcnological success of the rainbow- 
ladder kernel, manifest in the results of Tabled is unsurprising and, indeed, was 
to be expected. 



3.2 Comparison with Lattice Simulations 

The solution of the gap equation has long been of interest in grappling with 
DCSB in QCD and hence, in Figs. I depict the scalar functions characterising 
the renormalised dressed-quark propagator: the wave function renormalisation, 
Z(p 2 ), and mass function, M(p 2 ), obtained by solving Eq. (|15fl using Eq. 1)55(1 . 
The infrared suppression of Z(p 2 ) and enhancement of M(p 2 ) are longstanding 
predictions of DSE studies pfl. Indeed, this property of asymptotically free theo- 
ries was elucidated in Refs. JS| an d could be anticipated from studies of strong 
coupling QED |27j . The prediction has recently been confirmed in numerical 
simulations of quenched lattice-QCD, as is evident in the figures. 

It is not yet possible to reliably determine the behaviour of lattice Schwinger 
functions for current-quark masses that are a realistic approximation to those of 
the u- and d-quarks. A veracious lattice estimate of m n , f n , (qq) is therefore 
absent. To obtain such an estimate, Ref. [7j used the rainbow kernel described 
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herein and varied (D, ui) in order to reproduce the quenched lattice-QCD data. 
A best fit was obtained with 

D = (0.74 GeV) 2 , to = 0.3 GeV , (62) 

at a current-quark mass of 0.6 « 14 m u [Eq. 1)58(1 ] chosen to coincide with 

that employed in the lattice simulation. Constructing and solving the homoge- 
neous BSE for a pion-like bound state composed of quarks with this current-mass 
yields 

TO m„~14 m „ = Q_ 48 ym,~14 m „ = Q Q04 GeV . (63) 

The parameters in Eq. 1)62)1 give chiral limit results 0: 

/° = 0.068 GeV, -(gg)? GcV = (0.19 GeV) 3 , (64) 

whereas Eqs. lf57))l. (f5T)l give /° = 0.088 GeV. These results have been confirmed 
in a more detailed analysis [2 8) and this correspondence suggests that chiral and 
physical pion observables are materially underestimated in the quenched theory: 
I (qq) I by a factor of two and f v by 30%. 

The rainbow-ladder kernel has also been employed in an analysis of a tra- 
jectory of fictitious pseudoscalar mesons, all composed of equally massive con- 
stituents [2H| (The only physical state on this trajectory is the pion.) The DSE 
study predicts [301 

I^H-=2.2, (65) 

m H m=ms 

in agreement with a result of recent quenched lattice simulations The DSE 
study provides an intuitive understanding of this result, showing that it owes 
itself to a large value of the in-hadron condensates for light-quark mesons; e.g., 
(qq)i GeV = (—0.32 GeV) 3 0, and thereby confirms the large-magnitude con- 
densate version of chiral perturbation theory, an observation also supported by 
Eq. References [29132) also provide vector meson trajectories. 



3.3 Ab Initio Calculation of Meson Properties 

The renormalisation-group- improved rainbow- ladder kernel defined with Eq. 1)55)) 
has been employed to predict a wide range of meson observables, and this is re- 
viewed in Ref. @|. These results; e.g., those for vector mesons in Tabled are true 
predictions, in the sense that the model's mass-scale was fixed, as described in 
connection with Eq. I|57l) , and every element in each calculation was completely 
determined by, and calculated from, that kernel. 

A particular success was the calculation of the electromagnetic pion form 
factor, which is described in Refs. [33134) . The result is depicted in Fig.HJl wherein 
it is compared with the most recent experimental data £ 37 . It is noteworthy that 
all other pre-existing calculations are uniformly two - four standard deviations 
below that Q 2 F 7r (Q 2 ) data. 6 

6 The nature and meaning of vector dominance is discussed in Sec. 2.3.1 of Ref. |39j . 
Sec. 2.3 of Ref. |5] and Sec. 4.3 of Ref. [I]: the low-q 2 behaviour of the pion form factor 
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Fig. 6. Impulse approximation DSE prediction for q 2 F^(q 2 ) obtained in a parameter- 
free application of the renormalisation-group-improved rainbow-ladder truncation, Eqs. 
E3, E3- The data are from Refs. |3Sl36l3 l T| . (Adapted from Ref. |55).1 

In this connection one should also note that it is a model independent DSE 
prediction 0U] that electromagnetic elastic meson form factors display 



with calculable (lnq 2 /AQ CD ) d corrections, where d is an anomalous dimension. 
This agrees with earlier perturbative QCD analyses [41142] . However, to obtain 
this result in covariant gauges it is crucial to retain the pseudovector compo- 
nents of the Bethe-Salpeter amplitude in Eq. (J3SJl: F n , G^. (NB. The quark- 
level Goldberger-Treiman relations, Eqs. (|42|l . prove them to be nonzero.) With- 
out these amplitudes q 2 F(q 2 ) oc 1/q 2 . The calculation of Ref. 0U] sug- 
gests that the perturbative behaviour of Eq. ffify is unambiguously evident for 
q 2 > 15 GeV 2 . Owing to challenges in the numerical analysis, the ab initio calcu- 
lations of Ref. cannot yet make a prediction for the onset of the perturbative 
domain but progress in remedying that is being made • 

Another very instructive success is the study of ir-ir scattering, wherein 
a range of new challenges arise whose quiddity and natural resolution via a 
symmetry-preserving truncation of the DSEs is explained in Sec. 4.6 of Ref. 0], 
which reviews the seminal work of Refs. |44j . It is worth remarking, too, that 
with a systematic and nonperturbative DSE truncation scheme, all consequences 
of the Abelian anomaly and Wess-Zumino term are obtained exactly, without 



3.4 Heavier Mesons 

The meson spectrum contains [25^ four little-studied axial-vector mesons com- 
posed of u- and <i-quarks. They appear as isospin I = 0,1 partners (in the 

is necessarily dominated by the lowest mass resonance in the J PC = 1 channel. 
Any realistic calculation will predict that and also a deviation from dominance by 
the p-meson pole alone as spacelike-q 2 increases. 



q 2 F(q 2 ) = constant, q 2 > A% 



QCD j 



(66) 



fine tuning |:ffll4()l4KI46l4VI4s| . 
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manner of the w and p): Tlx (1170), & 1 (1235); and /i(1285), oi(1260), and differ 
in their charge-parity: J PC = 1 H for hi, bi, and J PC = 1 ++ for fi, a%. In 
the qq constituent quark model the bi is represented as a constituent-quark and 
-antiquark with total spin S = and angular momentum L = 1, while in the ai 
the quark and antiquark have S = 1 and L = 1. It is therefore apparent that 
in this model the bi is an orbital excitation of the tt, and the ai is an orbital 
excitation and axial- vector partner of the p. In QCD the J PC characteristics of a 
quark-antiquark bound state are manifest in the structure of its Bethe-Salpeter 
amplitude |1M| . This amplitude is a valuable intuitive guide and, in cases where 
a qq constituent quark model analogue exists, it incorporates and extends the 
information present in that analogue's quantum mechanical wave function. 

Three of the axial-vector mesons decay predominantly into two-body final 
states containing a vector meson and a pion: hi — > pit] bi — > lutt; a\ — ► p-n. 
With a J = 1 meson in both the initial and final state these three decays 
proceed via two partial waves (S, D), and therefore probe aspects of hadron 
structure inaccessible in simpler processes involving only spinless mesons in the 
final state, such as p — > nn. For example and of importance, in constituent- 
quark-like models the D/S amplitude ratio is very sensitive |49j to the nature 
of the phenomenological long-range confining interaction. 

The additional insight and model constraints that such processes can pro- 
vide is particularly important now as a systematic search and classification of 
"exotic" states in the light meson sector becomes feasible experimentally. I note 
that a meson is labelled "exotic" if it is characterised by a value of J PC which 
is unobtainable in the qq constituent quark model; e.g., the experimentally ob- 
served [nOl 7Ti(1600), a 1.6 GeV J PC = 1~ + state. Such unusual charge parity 
states are a necessary feature of a field theoretical description of quark-antiquark 
bound states |13j with BSE studies typically yielding 51 masses approximately 
twice as large as that of the natural charge parity partner and, in particular, a 
J PC = I" + meson with a mass ~ 1.5 GeV |52| . 

In appreciation of these points, Ref. [HSj used the simple DSE-based model 
of Ref. [SJ in a simultaneous study of axial-vector meson decays, p — > tttt de- 
cay, and the electroweak decay constants of the mesons involved. The results are 
instructive. It was found that the rainbow-ladder truncation is capable of simul- 
taneously providing a good description of these observables but that the D/S 
partial-wave ratio in the decays of axial-vector mesons is indeed very sensitive 
to details of the long-range part of a model interaction; i.e., to the expression 
of light-quark confinement. This is perhaps unsurprising given that the mass of 
each axial- vector meson mass is significantly greater than 2M (0); namely, twice 
the constituent-quark mass-scale. Unfortunately, more sophisticated calculations 
are lacking. This collection of experimentally well-understood mesons has many 
lessons to teach and should no longer be ignored. 
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Fig. 7. Quark mass function obtained as a solution of Eq. I|67|l using the rainbow 
truncation, discussed in connection with Eqs. 1121 - I15L and the interaction of Eq. 1551 
with current-quark masses, fixed at £ = 19 GeV: m Ut d{0 — 3.7 MeV, m 3 (£) = 82MeV, 
m c (C) = 0.58 GeV and m b (() = 3.8 GeV. The indicated solutions of M 2 (p 2 ) = p 2 define 
a Euclidean constituent-quark mass, Mf, which takes the values: = 0.56 GeV, 
M S E = 0.70 GeV, M C B = 1.3 GeV, M b B = 4.6 GeV. 

4 Heavy Quarks 

4.1 Features of the Mass Function 

The DSE methods described hitherto have been applied to mesons involving 
heavy-quarks 11 54 53j and in this case there is a natural simplification. To 
begin, one focuses on the fact that mesons, whether heavy or light, are bound 
states of a dressed-quark and -antiquark, with the dressing described by the gap 
equation, Eq. ft), written explicitly again here with the addition of flavour label, 
/(= u,d,s,c, b): 



Sfip)- 1 = i 1 -pA f (p 2 ) + B f (p 2 ) = A f (p 2 ) [ij-p + Mfip 2 )] (67) 




The other elements of Eq. (|68fl will already be familiar. 

The qualitative features of the gap equation's solution are known and typical 
mass functions, M*(p 2 ), are depicted in Fig.0 There is some quantitative model- 
dependence in the momentum-evolution of the mass-function into the infrared. 
However, with any Ansatz for the effective interaction that provides an accurate 
description of f Wt K and m^if, one obtains solutions with profiles like those 
illustrated in the figure. Owing to Eq. (fH5|) the ultraviolet behaviour is naturally 
fixed, namely, it is given by Eq. © for massive quarks and by Eq. lj4*3|l in the 
chiral limit. 
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It is apparent in the figure that as p 2 decreases the chiral-limit and u, d-quark 
mass functions evolve to coincidence. This feature signals a transition from the 
perturbative to the nonperturbative domain. Furthermore, since the chiral limit 
mass- function is nonzero only because of the nonperturbative DCSB mechanism, 
whereas the u, (i-quark mass function is purely perturbative at p 2 > 20GeV 2 , 
it also indicates clearly that the DCSB mechanism has a significant impact on 
the propagation characteristics of u, d, s-quarks. However, it is conspicuous in 
Fig. that this is not the case for the 6-quark. Its large current-quark mass al- 
most entirely suppresses momentum-dependent dressing so that Mf,(p 2 ) is nearly 
constant on a substantial domain. The same is true to a lesser extent for the 
c- quark. 

The quantity £/ := Mf /m/(£) provides a single quantitative measure of 
the importance of the DCSB mechanism; i.e., nonperturbative effects, in modi- 
fying the propagation characteristics of a given quark flavour. In this particular 
illustration it takes the values 



C 



f 



u, d 



150 10 2.2 1.2 



(69) 



which are representative: for light-quarks C q=u> d,s ~ 10-100; while for heavy- 
quarks C,Q =c ,b ~ I- They also highlight the existence of a mass-scale, M x , char- 
acteristic of DCSB: the propagation characteristics of a flavour with m,f(() < M x 
are significantly altered by the DCSB mechanism, while momentum-dependent 
dressing is almost irrelevant for flavours with m/(£) 3> M x . It is evident and 
unsurprising that M x ~ 0.2GeV~ Aqqd- Consequently one anticipates that 
the propagation of c, 6-quarks should be described well by replacing their mass- 
functions with a constant; viz., writing jll) 

Sq(p) = ^-rr:Q = c,6, (70) 

Z7 • p + Mq 

where Mq is a constituent-heavy-quark mass parameter. 7 

When considering a meson with an heavy-quark constituent one can proceed 
further, as in heavy-quark effective theory (HQET) [SS], allow the heaviest quark 
to carry all the heavy-meson's momentum: =: ra^v^ —: (Mq + Eh)vu,, and 
write 

i»zi!^ + o(JSL.*V (7i, 

V 2 k-v-E H {Mq M q ) 

where k is the momentum of the lighter constituent. It is apparent from the 
study of light-meson properties that in the calculation of observables the meson's 
Bethe-Salpeter amplitude will limit the range of |fe| so that Eq. I|71|l will only be 
a good approximation if both the momentum-space width of the amplitude, ujh, 
and the binding energy, Eh, are significantly less than Mq. 

In Ref. the propagation of c- and 6-quarks was described by Eq. 1)71)1 . 
with a goal of exploring the fidelity of this idealisation, and it was found to allow 

7 Although not illustrated explicitly, when Mf(p 2 ) w const., Af(p 2 ) ~ 1 in Eq. 16 711 . 
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for a uniformly good description of i?/-meson leptonic and semileptonic decays 
with heavy- and light-pseudoscalar final states. In that study, ujb s ~ 1.3 GeV 
and Eb } ~ 0.70 GeV, both of which are small compared with M& « 4.6 GeV in 
Fig. [7| Hence the accuracy of the approximation could be forseen. It is reasonable 
to suppose that ujd ~ ^>b and Ed »s Eb, since they must be identical in the limit 
of exact heavy-quark symmetry. Thus in processes involving the weak decay of a 
c-quark (M c ~ 1.3 GeV) where a Dj-meson is the heaviest participant, Eq. Q71|l 
must be inadequate; an expectation verified in Ref. |55|. 

The failure of Eq. I|71(l for the c-quark complicates or precludes the develop- 
ment of a common understanding of Dt- and £>/-meson observables using such 
contemporary theoretical tools as HQET and light cone sum rules. However, as 
shown in Ref. ^I] an d I will illustrate, the constituent-like dressed-heavy-quark 
propagator of Eq. (|70l can still be used to effect a unified, accurate simplification 
in the study of these observables. 

4.2 Leptonic Decays 

Pseudoscalar Mesons. The leptonic decay of a pseudoscalar meson, P(p), is 
described by the matrix element (Sec. I2.4JI 

r A 

fpPv ■= <0|Q(T p ) t 7ai75 Q\P(p)) = trZ 2 J {T p ) T lbl , XP {k-p), (72) 

where Q = column(u, d, s, c, b) and here I have adopted a charged particle nor- 
malisation, which yields results for fp a factor of y/2 larger than Eq. l(5^|) and 
is conventional in studying heavy-quark systems. 

In Eq. Xp is the meson's Bethe-Salpeter wave function, related to its 

amplitude, Jp, via Eq. (|3 1|) and normalised canonically as described in connec- 
tion with Eq. <|30|l . Using Eq. I|71|l . it follows from the canonical normalisation 
condition that 

Gp{k;p) := —. - — r P (k; p) < oo, m P -> oo; (73) 

i.e., Qp{k;p) so-defined is mass-independent in the heavy-quark limit. Using this 
result plus Eq. J7TJ one finds from Eq. (5^ 

fp °c —j- — , mp-*oo. (74) 
V m P 

Equation (|74|) is a model-independent result and a well-known general con- 
sequence of heavy-quark symmetry |56j . However, the value of the hadron mass 
at which this behaviour becomes evident is unknown. It is clear from Table [3 
that, experimentally, 



U = 131 MeV < f K = 160 MeV 



(75) 
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Furthermore, direct DSE studies following the method described in Sec. l3.3l show 
that for pseudoscalar mesons uf, composed of a single u, d-quark and an anti- 
quark of mass m/, fp(mp) is a monotonically increasing concave-down function 
on mp € [0, 0.9] GeV, where mp is the calculated mass of this composite system, 
and likely on a larger domain |32| . On the other hand, numerical simulations of 
quenched lattice-QCD indicate [57] 

f D = 200±30MeV > f B = 170±35MeV. (76) 

In simulations of lattice-QCD with two flavours of sea quarks both of these decay 
constants increase in magnitude but there is no sign that the ordering is reversed 
|57I58| . The information in Eqs. Ij75(l . Ij76(l is depicted in Fig. El This and analysis 
to be reviewed subsequently suggest that D-mesons lie outside the domain on 
which Eq. (|74|l is a reliable tool. 



Vector Mesons. The leptonic decay constant, fy, for a vector meson with 
mass My is given in Eq. (|61|l and adapting the analysis that leads to Eq. Q74|l 
one finds readily 

f v cx — , My > oo , (77) 

which again is a model-independent result. Moreover, since the pseudoscalar and 
vector meson Bethe-Salpeter amplitudes become identical in the heavy-quark 
limit, it follows that JT] 

fv = fp > My = mp , in the limit mp — > oo . (78) 



4.3 Heavy-Meson Masses 

More can be learnt from the pseudoscalar meson mass formula in Eq. I|52l) . Using 
Eq. (|74|) . and applying to Eq. (|54|l the analysis from which it follows, one obtains 

— (Ql)c = constant , as mp —> oo (79) 
and consequently |l(Jlllj 

mp oc rriQ , ttiq — > oo , (80) 

where ttiq is the renormalisation-group-invariant current-quark mass of the fia- 
vour-nonsinglet pseudoscalar meson's heaviest constituent. This is the result one 
would have guessed from constituent-quark models but here I have outlined a 
direct proof in QCD. 

Equation l|47|) is thus seen to be a single formula that unifies the masses of 
light- and heavy-quark mesons. This aspect has been quantitatively explored 
using the rainbow-ladder kernel described in Sec. El with the results illustrated 
in Fig. El Therein the calculated mass of a uf pseudoscalar meson is plotted 
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Fig. 9. Pseudoscalar uq meson's mass as a function of m q (Q, £ = 19 GeV, with a 
fixed value of m u (C,) corresponding to ra u (lGeV) = 5.5 MeV, Eq. 158H (solid line). 
The experimental data points are from Ref. |25| as are the errors assigned to the 
associated heavy-quark masses. A straight line is drawn through the K, D, B masses 
(dashed curve). (Adapted from Ref. See also Ref. |32|.1 



as a function of m/(C), with m u {Q fixed at the value in Eq. i|58[l. The DSE 
calculations are depicted by the solid curve, which is (2H| (in MeV) 



m P = 83 + 500V^ + 310 X, X = m c q /A Q cr>. 



(81) 
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The curvature appears slight in the figure but that is misleading: the nonlinear 
term in Eq. I jSljl accounts for almost all of (the Gell-Mann-Oakes-Renner 
relation is nearly exact for the pion) and 80% of mx- NB. The dashed line in 
Fig. El fits the K, D, B subset of the data exactly. It is drawn to illustrate how 
easily one can be misled. Without careful calculation one might infer from this 
apparent agreement that the large-m g limit of Eq. (|47|) is already manifest at 
the s-quark mass whereas, in reality, the linear term only becomes dominant for 
m q > 1 GeV, providing 50% of mp, and 67% of mp- The model predicts, via 
Eq. {52), m;l GeV = 1.1 GeV and mJ GcV = 4.2 GeV, values that are typical of 
Poincare covariant treatments. 



4.4 Semileptonic Transition Form Factors 

Pseudoscalar meson in the final state. The transition: P\{pi) — > ^2(^2) t-V , 
where P\ represents either a B- or D-meson and P2 can be a D, K or 7r, is 
described by the invariant amplitude 

A(P 1 ^P 2 iv) = ^ V n %(1 - 75)^ M^( Pl , P2 ) , (82) 

where Gf = 1.166 x 10~ 5 GeV -2 , Vfj is the relevant element of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix, and the hadronic current is 

M^( Pl ,p 2 ) := (P2(p2)\f'^f\Pi(pi)) = /+(*)(Pi +P2)»+f-(t)q», (83) 

with t := — q 2 = — (pi — P2) 2 ■ The transition form factors, f±(t), contain all 
the information about strong-interaction effects in these processes, and their 
accurate calculation is essential for a reliable determination of the CKM matrix 
elements from a measurement of the decay width (t± :— (m Pl ± rap 2 f): 

r(p 1 ^p 2 £v) = -^\v rf \ 2 -± r [*~ dt\f + (t)\ 2 [{t + -t)(t_-t)f 2 . (84) 
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The related study of light-meson initial states is described in Refs. |60|. 

Vector meson in the final state. The transition: P( P i) — > V\(p 2 )£v, with 
P either a B or D and V\ a D*, K* or p, is described by the invariant amplitude 

A(P -^ Vx i v ) = 9jL Vflf U 1im{ i _ l5> e A (p2) M PVx {pup2) f (85) 

in which the hadronic tensor involves four scalar functions 

4(P2) M^( Pl ,p 2 ) = €^ (m P + My) A 1 (t) + ( Pl +p 2 ), e x ■ q A f] 

p p nip + My 

A Mt) A 2V(t) 

+ 9 - £ ' 9 m P + My + ^ e « PlaP20 m P + M v - 

(86) 
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Introducing three helicity amplitudes 

H ± = (rn P + My) A 1 (t) T A *( ro ^' t ) V (t) , (87) 

m P + My 



where A(mp, My, t) = [i + — t) [t- —t],t±= (rap ± My) 2 , the transition rates 

TT = I 5^L-|V'„| 2 tAJ(^,Mj, f )| ff± ,„(,)r, (88) 

from which one obtains the transverse and longitudinal rates 

dr T dr+ dP_ /"*- dTT , x 



dt dt dt ' J Q dt 

dr L dr „ r*- „ dr L 



dt dt> rL = J Q d< ^F< < 91 > 

wherefrom the total width r = + Tl- The polarisation ratio and forward- 
backward asymmetry are 

a = 2^-l, A ¥B = \ r -^±. (92) 
4.5 Impulse Approximation 

As explained and illustrated in Ref. the impulse approximation is accurate 
for three point functions and applied to these transition form factors it yields 



2N c tx D / f x (k;-p2) S g (k 2 )iOf(k 2> k 1 ) Sq(ki) r P (k;pi) S q >(k), (93) 

Jk 

wherein the flavour structure is made explicit, ki t2 = k +Pi,2 and: 

H* x=p >fa,P2) = M^( Pl , P2 ) , (94) 
K x - v \ PuP2 ) = et( P2 )M^( Pl , P2 ) ; (95) 

r x= y*{k; P ) — € X ( P ) ■ r v (k; P ); and OjQ(k 2 ,ki) is the dressed-quark-W-boson 
vertex, which in weak decays of heavy-quarks is well approximated by 54 55 

Of(k 2 ,k 1 )= 1 ,(l- l5 ) (96) 

because Aq( P 2 ) w const, and Mq(p 2 ) w const, for heavy-quarks (recall Fig. H) 
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Quark Propagators. It is plain that to evaluate H^ lX (pi,P2) a specific form 
for the dressed-quark propagators is required. Equation l|70[) provides a good 
approximation for the heavier quarks, Q — c,b, as explained in Sec. 14. II and this 
was used in Ref. with Mq treated as free parameters. 
For the light-quark propagators: 



S f(?) = -ij-pVv(P 2 )+ a s(P 2 ) 



1 



i~fpA f (p 2 ) + Bf(p 2 ) 



(97) 



Ref. JT] assumed isospin symmetry and employed the algebraic forms intro- 
duced in Ref. [211] , which efficiently characterise the essential features of the gap 
equation's solutions: 



a f s {x) = 2m f T(2(x 



fn 2 )) 



a v (x) 



2{x- 



1 + e^ 2 ^" 1 /) 



b{T{ex) 



2(x- 



m 2 f ) 2 



F(y) = (1 — e v )/y, x = p 2 /X 2 : nrif = mj/X; and (Jg{x) 



-J 



(98) 
(99) 

(x) = 

X 2 <jy(p 2 ), with A a mass scale. The parameters are the current-quark mass, in, 
and 6o,i,2,3; about which I shall subsequently explain more. 

This algebraic form combines the effects of confinement 8 and DCSB with 
free-particle behaviour at large, spacelike p 2 . One characteristic of DCSB is the 
appearance of a nonzero vacuum quark condensate and using this parametrisa- 
tion in Eq. (|5U)l yields 



Xa f s (p 2 ) 



= X 



4?r 2 b^b^ 



In 



c 2 



A 2 

QCD 



(100) 



The simplicity of this result emphasises the utility of an algebraic form for the 
dressed-quark propagator. That utility is amplified in the calculation of a form 
factor, which requires the repeated evaluation of a multidimensional integral 
whose integrand is a complex- valued function, and a functional of the propagator 
and the Bcthc-Salpctcr amplitudes. 



Bethe-Salpeter Amplitudes. An algebraic parametrisation of the Bethe- 
Salpeter amplitudes also helps and the quark-level Goldberger-Trciman relation, 
Eq. (|42p. suggests a form for light pseudoscalar mesons: 

r P (k;p) =i l5 £ P (k 2 ) = i l5 -L B P {k 2 ) , P = ir,K, (101) 
fp 

8 The representation of S(p) as an entire function is motivated by the algebraic solu- 
tions of Eq. in Refs. The concomitant violation of the axiom of reflection 
positivity is a sufficient condition for confinement, as reviewed in Sec. 6.2 of Ref. 0, 
Sec. 2.2 of Ref. and Sec. 2.4 of Ref. 0. 
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where Bp := B u \ bU ^ bP , obtained from Eq. Q97|). and fp = fp/^/2 because in 

this Section I use the / OT = 131 MeV normalisation. b^' K are two additional pa- 
rameters. This Ansatz omits the pseudovector components of the amplitude but 
that is not a material defect in applications involving small to intermediate mo- 
mentum transfers 0U| . Equations (|52ll , 154|) , (|101|) yield the following expression 
for the 7r- and K -meson masses: 

fp nip = -(m u + rrifp) (qq)i GcV 2 , (102) 

where — ma, rrifK — m s , and 

In studies of the type reviewed in Sec. this in-hadron condensate takes values 

(99>iGoV 2 ~ l-05<Uu) 1G eV 2 and (<?<?}f G cV 2 ~ i'S^lGcV 2 - 

Employing algebraic parametrisations of the light vector meson Bethe-Salpeter 
amplitudes is also a useful expedient and that approach was adopted in Rcfs. 
|1 1154155] . Indeed, sophisticated calculations of light vector meson properties 
based on the rainbow-ladder truncation did not exist at the time of those stud- 
ies, although it was clear that a given vector meson is narrower in momentum 
space than its pseudoscalar partner, and that for both vector and pseudoscalar 
mesons this width increases with the total current-mass of the constituents. 
These qualitative features were important in the explanation of meson electro- 
production cross sections and electromagnetic form factors [SH], and can be 
realised in the simple expression 

r /T(*5P) = jjv +Pm J^t) V(k 2 ) , (104) 

where f(k 2 ) = 1/ (l+fc 4 /wy) with uiy a parameter and Af v fixed by the canonical 
normalisation condition. One expects: wjf. ps 1.6 to p |63) . 

In connection with the impulse approximation to scmilcptonic transition form 
factors it remains only to fix the heavy- meson Bethe-Salpeter amplitudes. In this 
case, too, algebraic parametrisations offer a simple, attractive and expeditious 
means of proceeding and that again was the approach adopted in Ref. |11| . 
Therein heavy vector mesons were described by Eq. (|104fl . with <p(k 2 ) — > ipij(k 2 ), 
and heavy pseudoscalar mesons by its analogue: 

r P (k-p) = ^pi ToVH {k 2 ) 1 (105) 

where (pij(k 2 ) = exp (— k 2 /^ 2 ^) . The amplitudes are again normalised canoni- 
cally. Such a parametrisation naturally introduces additional parameters; viz., 
the widths. The number is kept at two by acknowledging that Bethe-Salpeter 
amplitudes for truly heavy-mesons must be spin- and flavour-independent and 
assuming therefore that lob = ^b* = wb, and ujp> = lub>* = ujd 3 - 
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4.6 Additional Decay Processes 

Many more decays were considered in Ref. , with the goal being to determine 
whether a unified description of light- and heavy-meson observables is possible 
based simply on the key DSE features of quark dressing and sensible bound 
state amplitudes. For example, there are experimental constraints on radiative 
decays H* — > H~/, where H = D/ s \, Br s \, and so these widths, Fh*^Hi\ were 
calculated. The strong decays H* — > H tt were also studied. They can be charac- 
terised by a coupling constant gH*Hir, which is calculable even if the process is 
kincmatically forbidden, as is B* — > Bit. Lastly, the width for the rare flavour- 
changing neutral current process B — ► K*j, which proceeds predominantly via 
the local magnetic penguin operator j^l] and can be characterised by a coupling 
9bk*~/, was calculated because data exists and this process might be expected 
to severely test the framework since it completely exceeds the scope of previous 
applications. 



4.7 Heavy-Quark Symmetry Limits 

Equation l|71|) . and Eq. Ij73() and its natural analogues, can be used to elucidate 
the heavy-quark symmetry limit of the impulse approximation to any process 
and many were made explicit in Refs. [11154155] . I will only recapitulate on the 
most straightforward three-point case; namely, the semileptonic heavy — > heavy 
transitions. From Eqs. (|83|l . (|93|l one obtains 



f±(t)=T ± t; f (i 



where: W 



f 



4tt 2 



mp 2 ±mp 1 
2^/mp 2 m Pl 

1 



dr- 



ii 



duip H (z w ) 







(106) 



a f s (z w ) + J ^-<j f v {z w ) 



ran n 



f 



1 + 2t(1-t)(w-1), z w 

i n c r 

m H 47t 2 J 



■ N 2 P 



--u-2E H ^MW\ 
du(p 2 H (z) {a f s (z) + ^ua{,(z)j 



(107) 



with z = u — 2Eh\/u, f labelling the meson's lighter quark and all dimensioned 
quantities expressed in units of the mass-scale, A; and 



w 



l P2 



t 



-v Pl ■ Vp 2 



(108) 



2mp 1 mp 2 

The canonical normalisation of the Bethe-Salpeter amplitude automatically en- 
sures 

£ f (w = 1) = 1 (109) 

and from Eq. (|106|l follows 



P 



dw 



> 



1 



(110) 
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Semileptonic transitions with heavy vector mesons in the final state, de- 
scribed by Eqs. (|86l) and (|93() . can be analysed in the same way, and that yields 

A^t) = -L i (i + w) £ /H , A 2 (t) = -A 3 (t) = V(t) = T+ Zfiw) . (Ill) 

Equations (jl06[l . (|1 1 lfl are exemplars of a general result that in the heavy- 
quark symmetry limit the semileptonic Hf — > H'^ transitions are described by a 
single, universal function: £/(u>) |65|. In this limit the functions 

R l {w) := (1 - t/t + ) Z|L , R 2 (w) := (1 - t/t + ) M± (112) 

are constant (=1), independent of w. 



4.8 Survey of Results for Light- and Heavy-Meson Observables 

With every necessary element defined, the calculation of observables is a straight- 
forward numerical exercise. The algebraic Ansatze described above involve ten 
parameters plus four current-quark masses and in Ref. |11| they were fixed via a 
X 2 -fit to the iVobs = 42 heavy- and light-meson observables in Table |3| a process 
which yielded [HEj 





rrif b{ b^ 2 




. ,GeV 


. ,GcV 


M£eV 


u 
s 


0.00948 2.94 0.733 tt 
0.210 3.18 0.858 K 


0.204 p 
0.319 K* 


0.515 D 
0.817 B 


1.81 c 
1.81 b 


1.32 (113) 
4.65 



with x 2 /d.o.f = 1.75 and xV-^obs = 1-17. The dimensionless u,s current-quark 
masses correspond to m u = 5.4 MeV, m s — 119 MeV, and M® — 0.36 GeV, 
Mf = 0.49 GeV. Furthermore, lok*I^p = 1-59, which is identical to the value 
in Ref. gig. 

It is evident that the fitted heavy-quark masses are consistent with the esti- 
mates in Ref. |25| and hence that the heavy-meson binding energy is large: 

E D := m D - M c = 0.67 GeV, E B := m B - M b = 0.70 GeV. (114) 

These values yield Ed/M c = 0.51 and Eb/Mi, = 0.15, which furnishes another 
indication that while an heavy-quark expansion is accurate for the 6-quark it 
will provide a poor approximation for the c-quark. This is emphasised by the 
value of u>d = lob, which means that the Compton wavelength of the c-quark is 
greater than the length-scale characterising the bound state's extent. 

With the parameters fixed, in Ref. values for a wide range of other 
observables were calculated with the vast majority of the results being true 
predictions. The breadth of application is illustrated in Table and in Fig. 1101 
which depicts the calculated i-dependence of B — > 7r , p semileptonic transition 
form factors. I note that now there is a first experimental result for the D* + 
width [72j: r D .+ = (96±4±22)keV, g D . Dv = 17.9 ±0.3 ±1.9. Its confirmation 
and the gathering of additional information on c-quark mesons is crucial to 
improving our knowledge of the evolution from the light- to the heavy-quark 
domain, a transition whose true understanding will significantly enhance our 
grasp of nonperturbative dynamics. 
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Table 3. The 16 dimension-GeV (upper panel) and 26 dimensionless (lower panel) 
quantities used in the X 2_n t of Ref. The weighting error was the experimen- 

tal error or 10% of the experimental value, if that is greater, which accounts for a 
realistic expectation of the model's accuracy. The light-meson electromagnetic form 
factors were calculated in impulse approximation 39 40 167) and £(u>) was obtained 
from f+^ D (t) via Eq. (TUTSI . The values in the "Obs." column were taken from Refs. 
l6l^l36K7l68l69TrWTl . (Adapted from Ref. [TT].) 





Obs. 


Calc. 




Obs. 


Calc. 


A 


0.131 


0.146 


m w 


0.138 


0.130 


Ik 


0.160 


0.178 


m K 


0.496 


0.449 


(uu) 1/z 


0.241 


0.220 


<-> 1/3 


0.227 


0.199 




0.245 


0.255 




0.287 


0.296 


f P 


0.216 


0.163 




0.244 


0.253 




0.151 


0.118 


•Tk * (Kir) 


0.051 


0.052 


fo 


0.200 ± 0.030 


0.213 


fo, 


0.251 ± 0.030 


0.234 


fs 


o i 7d + n rn^ 

u.ilu zn u.uoj 


n 1 89 




9 rn + n f>9 


9 8fi 


fr°(o) 


0.73 


0.58 


fnr n 


0.44 ± 0.004 


0.44 


^(S.SGcV 2 ) 


0.097 ± 0.019 


0.077 


B(B -> D*) 


0.0453 ± 0.0032 


0.052 


p 2 


1.53 ± 0.36 


1.84 


B->D* 

a 


1.25 ± 0.26 


0.94 


£(l-l) 


0.86 ± 0.03 


0.84 




0.19 ± 0.031 


0.24 


£(1.2) 


0.75 ± 0.05 


0.72 


B(B -> tt) 


(1.8±0.6) xlo -4 


2.2 


£(1-3) 


0.66 ± 0.06 


0.63 


rB^rr 

J + (14.9 GeV 2 ) 


0.82 ± 0.17 


0.82 


£(1-4) 


0.59 ± 0.07 


0.56 


fB — '7T 

■'+ (17.9 GeV 2 ) 


1.19 ± 0.28 


1.00 


£(1-5) 


0.53 ± 0.08 


0.50 


f B — '7T 

J+ (20.9 GeV 2 ) 


1.89 ± 0.53 


1.28 


B(B -> D) 


0.020 ± 0.007 


0.013 


B(B - p) 


(2.5±0.9) xl0 -4 


4.8 


B(D -f K*) 


0.047 ± 0.004 


0.049 




0.73 


0.61 


V(0) D - K * 

^i(O) 

r T 

a 2 (q) d - k * 
Ai(0) 


1.89 ± 0.25 
1.23 ± 0.13 


1.74 
1.17 




0.73 

23.0 ± 5.0 


0.67 
23.2 


0.73 ± 0.15 


0.87 




10.0 ± 1.3 


11.0 



Fidelity of heavy-quark symmetry. The universal function characterising 
semileptonic transitions in the heavy-quark symmetry limit, £(uj) introduced in 
Sec. 14.71 can be estimated most reliably from B — > D, D* transitions. Using Eq. 
1)106)1 to infer this function from f+^ D (t), one obtains 



£'+(!) = 1.08, 



(115) 
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Table 4. Predictions in Ref. for a selection of observables. The "Obs." values 
are extracted from Refs. ,25 57 68 72 73 . t max is the maximum momentum transfer 
available in the process identified and cj ma x = ii>(t maj[ ) calculated from Eq. I|108(l . 
(Adapted from Ref. [Hp 



fKr K 

Q ' p'KTl 

9p 

fo* (GeV) 
fo; (GeV) 
fB, (GeV) 
/a. (GeV) 
fa* (GeV) 
fo./fo 

f B Jf B 
fD'/fB 

fa*/ fa 
R?-* D * (uw) 



B(D+ - 
B(D° ~ 

A^ K '{Q: 

rD—*K* 



(0) 



v u "* ■ (0) 

B(D° -» tt) 



A') 



BCD 

Af^*(l) 
<^(1) 



Obs. 

0.472 ± 0.038 
6.05 ± 0.02 
6.41 ± 0.06 
5.03 ± 0.012 



0.195 ± 0.035 



1.10 ± 0.06 
1.14 ± 0.08 



1.30 ± 0.39 



0.037 ± 0.002 

0.56 ± 0.04 
0.39 ± 0.08 

1.1 ± 0.2 

0.103 ± 0.039 

1.2 ± 0.3 



Calc. 

0.46 

5.27 

5.96 

5.27 

0.290 

0.298 

0.194 

0.200 

0.209 

1.10 

1.07 

1.36 

1.10 

1.32 

1.23 

0.032 

0.036 

0.46 
0.40 
0.80 
0.098 

1.10 

1.72 
2.08 



r D , (MeV) 
r D , + (keV) 
r DtDsl (MeV) 
r B , +B+1 (keV) 
r B ,o BOj (keV) 
r B1Bsl (keV) 
B(£>*+-> D+ir°) 
B(D*+^ D°ir+) 
B(D*+^ D+7) 
B(D*°^D°tt°) 
B(D* ^ D°j) 
B(B -> A* 7 ) 



(1) 



R"^ U (Wmax) 



B(D 



B(D -> A*) 

"1 V^max J 



2 V C max ) 



A 



v V max 



+ V'rnax ) 



rD—>ir fiD — >7T\ 
J+ v^max ) 

R\ (ifimai) 
i?f^ P (TOmax) 



Obs. 

(0.19 ± 0.05) 2 

< 2.1 

96 ± 4 ± 22 

< 1.9 



0.306 ± 0.025 
0.683 ± 0.014 
011 + 021 

U.UXi. _ .007 

0.619 ± 0.029 
0.381 ± 0.029 
(5.7 ± 3.3) 10 - E 
0.64 ± 0.29 



0.044 ± 0.034 

0.66 ± 0.05 
0.46 ± 0.09 
1.4 ± 0.3 

1.31 ± 0.04 



Calc. 

(0.10) 2 

0.020 

37.9 

0.001 

0.030 

0.015 

0.011 

0.316 

0.683 

0.001 

0.826 

0.174 

11.4 

1.04 

0.98 

1.03 

0.065 

0.47 
0.44 
0.92 

1.11 
2.18 

1.74 
2.03 



which is a measurable deviation from Eq. Q109|l . The ratio (w)jij + (0) is 
depicted in Fig. 1111 wherein it is compared with two experimental fits |69| : 

£{w) = l-p 2 (w-l), p 2 = 0.91 ±0.15 ±0.16, (116) 
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Fig. 10. Semileptonic B — > n and B — > p form factors with, for comparison, data 
from a lattice simulation |70| and a vector dominance, monopole model: f+^ n (t) — 
0.46/(1 - t/m 2 B ,), m B * = 5.325 GeV {short-dashed line). (Adapted from Ref. |H).j 



exp 



(l-2p 2 



1 



p 2 = 1.53 ±0.36 ±0.14. (117) 



The evident agreement was possible because Ref. QT] did not employ the heavy- 
quark expansion of Eq. Ij71|l . in particular and especially not for the c-quark. 
The calculated result (solid curve) in Fig. 1111 is well approximated by 



1 



= 1.98. 



(118) 



Equations l|lll|) were also used in Ref. JT] to extract £(w) from B — > D*. 
This gave £ Al (l) = 0.987, £ A2 (1) = 1.03, £ v (l) = 1.30, an ^-dependence well- 
described by Eq. (|TT5|) but with p 2 Ai = 1.79, p\ = 1.99, p% = 2.02, and the 
ratios, Eqs. dU), Ri(l) / Ri(w max ) = 1.08, R 2 (l) / R 2 (w max ) = 1.06. 

This collection of results indicates the degree to which heavy-quark symmetry 
is respected in b — > c processes. Combining them it is clear that even in this case, 
which is the nearest contemporary realisation of the heavy-quark symmetry limit, 
corrections of < 30% must be expected. In c — * s , d transitions the corrections 
can be as large as a factor of two, as evident in Tabled 



5 Epilogue 

This contribution provides a perspective on the modern application of Dyson- 
Schwinger equations (DSEs) to light- and heavy-meson properties. The keystone 
of this approach's success is an appreciation and expression of the momentum- 
dependence of dressed-parton propagators at infrared length-scales. That depen- 
dence is responsible for the magnitude of constituent-quark and -gluon masses, 



Unifying aspects of light- and heavy-systems 



37 



1.00 

g 0.80 

| 

0.60 
0.40 




1.0 1.2 1.4 1.6 

w 

Fig. 11. Calculated form of £(w) (solid line) compared with experimental analyses: 
linear fit from Ref. |69| . Eq. 11161 (short-dashed line); nonlinear fit from Ref. |69|. 
Eq. I|117|l (long-dashed line). The two lighter dotted lines are the nonlinear fit of Ref. 
|69| evaluated with the extreme values of p 2 : upper line, p 2 = 1.17 and lower line, 
p 2 = 1.89. The data points are from Ref. |71]. (Adapted from Ref. [TT].1 

and the length-scale characterising confinement in bound states; and is now 
recognised as a fact. 

It has recently become clear that the simple rainbow-ladder DSE truncation 
is the first term in a systematic and nonperturbative scheme that preserves the 
Ward-Takahashi identities which express conservation laws at an hadronic level. 
This has enabled the proof of exact results in QCD, and explains why the trun- 
cation has been successful for light vector and flavour nonsinglet pseudoscalar 
mesons. Emulating more of these achievements with ab initio calculations of 
heavy-meson properties is a modern challenge. 

However, at present, the study of heavy-meson systems using DSE methods 
stands approximately at the point occupied by those of light-meson properties 
seven - eight years ago. A Poincare covariant treatment exploiting essential 
features, such as propagator dressing and sensible bound state Bethe-Salpeter 
amplitudes, has been shown capable of providing a unified and successful de- 
scription of light- and heavy-meson observables. The goal now is to make the 
case compelling by tying the separate elements together; namely, relating the 
propagators and Bethe-Salpeter amplitudes via a single kernel. I am confident 
this will be accomplished, and the DSEs become a quantitatively reliable and 
intuition building tool as much in the heavy-quark sector as they are for light- 
quark systems. 

While a more detailed understanding will be attained in pursuing this goal, 
certain qualitative results established already are unlikely to change. For ex- 
ample, it is plain that light- and heavy-mesons are essentially the same, they 
are simply bound states of dressed-quarks. Moreover, the magnitude of the b- 
quark's current-mass is large enough to sustain heavy-quark approximations for 
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its propagator and the amplitudes for bound states of which it is a constituent. 
In addition, and unfortunately in so far as practical constraints on the Standard 
Model are concerned, the current-mass of the c-quark is too small to validate an 
heavy-quark approximation. 
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